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Abstract. We review basic properties of advection-dominated accretion 
flows (ADAFs) and their applications to astrophysical systems ranging 
from Galactic binary systems to galactic nuclei. A new classification 
*n . scheme for low-luminosity, X-ray bright galactic nuclei is highlighted. 

^ | Some outstanding unresolved issues are discussed. 

in 
o 

o\ 

^ ' 1. Introduction 

Q-f It has recently been recognized that rotating accretion flows with low radia- 

q ■ five efficiency are applicable to a wide range of astrophysical systems including 

| Galactic X-ray transients and active galactic nuclei. In the recently discussed 

models for such sources, the radiative efficiency becomes low because the gravita- 
tional binding energy dissipated during infall of accreted matter is not efficiently 
radiated away due to low gas densities but stored within the flows and radially 
advected inward. Viscous torque transports angular momentum while generat- 
■ ing heat through dissipation. Radiative efficiency is basically determined by the 

fraction of the viscously dissipated energy which goes into radiation (Narayan 
& Yi 1995b). In the optically thin limit, which occurs when the density of the 
accretion flow falls below a certain limit, the radiative efficiency becomes small 
due to the long time scales for electron-ion energy exchange process and rele- 
vant cooling processes. The first discussions of such flows are found in Ichimaru 
(1977), Rees et al. (1982), and references therein. 

The accretion flows around compact objects could be classified into several 
types (e.g. Narayan et al. 1998b, Frank et al. 1992). (i) Geometrically thin 
disks radiate away a large fraction of the viscously dissipated energy and due 
to high density, the optical depth for outgoing radiation is large. The disk 
temperature is relatively low and hence the internal pressure support is small, 
which results in the small geometrical thickness (e.g. Frank et al. 1992 for 
review). In these flows, the cooling rate Q~ is balanced by the heating rate Q + , 
Q + = Q~ while the electron temperature T e equals the ion temperature Tj. (ii) 
Two-temperature, Shapiro-Lightman-Eardley (1976) type accretion flows also 
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maintain the energy balance as Q + = Q~ while T; t ~ T v ; ir ^> T e where T v ; ir is the 
usual virial temperature. These flows are thermally unstable (Piran 1978, cf. 
Rees et al. 1982). (hi) Slim disk or optically thick advection-dominated accretion 
flows (ADAFs) occur when accretion rates typically exceed the Eddington rate 
and the optical depth is large. In this limit, photons diffuse out on a time scale 
longer than the radial inflow time scale (Abramowicz et al. 1988 and references 
therein). These flows have Q + > Q~ and Tj = T e . (iv) Optically thin ADAFs 
(Ichimaru 1977, Rees et al. 1982, Narayan & Yi 1994, 1995b, Abramowicz et al. 
1995) are relevant for substantially sub-Eddington accretion rates. The optical 
depth in these flows are typically less than unity and the radiative cooling rate 
is much smaller than the heating rate, i.e. Q + > Q~ , while the ion temperature 
Tj ~ T v i r is much higher than the electron temperature T e . These flows are most 
interesting in systems which have relatively low luminosities and high emission 
temperatures (Narayan & Yi 1995b). 

In this review, we focus on the optically thin ADAFs and discuss their 
applications to various astrophysical systems. 



2. Basics of ADAFs: A Simple Version 
2.1. Basic Equations 

Following Narayan k, Yi (1994, 1995ab), we adopt the following notations for a 
steady, axisymmetric, rotating accretion flow; R = cylindrical radius from the 
central star, p = gas density, H = thickness or vertical scale height of the flow, 
v = radial velocity, f2 = angular velocity, tlx = (GM/R 3 ) 1 / 2 = Keplerian an- 
gular velocity, c s = isothermal sound speed, v = kinematic viscosity coefficient, 
T = is temperature of the gas, and s = specific entropy of the gas. Then, the 
basic conservation equations for mass, radial momentum, angular momentum, 
and energy become respectively, 

pRHv = constant, (1) 

„ rT d(nR 2 ) i / „« TT dn\ 

pvT—= q + -q-^fq + (4) 

where q + = v pR 2 '(dtt / 'dR) 2 = viscous dissipation rate per unit volume, q = 
radiative cooling rate per unit volume, and pvT(dsfdR) = q adv = radial ad- 
vection rate per unit volume. Therefore, the energy equation becomes simply 
q adv = q + — q~ = fq + which defines the advection fraction /. When advec- 
tion cooling per unit volume dominates, q~ q + , / ~ 1. Even when / differs 
substantially from unity, most of dynamical calculations assume a constant f 
rather than solving the full energy equation (or assumes a simple cooling such 
as bremsstrahlung) . Viscosity coefficient is specified by the a prescription (e.g. 
Frank et al. 1992) v = ac s H = ac 2 /Qx where a is a constant often assumed to 
be in the range 0.01 — 1, 
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The accretion flows are classified into three types according to the relative 
importance of the terms in the energy equation (Narayan et al. 1998). (i) q + ~ 
q~ 3> q adv : energy balance is maintained between viscous heating and radiative 
cooling, which corresponds to high-efficiency flows such as geometrically thin 
disks, (ii) q + ~ q adv 3> q~: radiative loss is negligible and luminosity is very 
low. (hi) g + « f « _qadv. v i scous heating is negligible and thermal energy 
of the flow is converted to radiation as in the cooling flows. Infall of matter is 
driven by pressure loss as gas cools. 

2.2. Self-Similar Solution 

The dynamical equations derived admit a self-similar solution for / =constant 
as shown in Narayan & Yi (1994) and Spruit et al. (1987); 
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where v K = R^K = (GM/R) 1 / 2 , e' = e/f, e = (5/3 - 7 )/( 7 - 1), and g(a,e') = 
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1. The specific heat ratio 7 = 4/3 — 5/3. 
For a = 0.01 - 0.3, a 2 < 1 and / w 1 gives 
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The self-similar solution reveals the basic properties of the ADAFs. (i) The 
radial accretion time scale is much shorter than that of the thin disk, (ii) Sub- 
Keplerian rotation occurs due to large internal pressure support, (iii) Vertical 
scale height H ~ c s /Qk ~ R- Moreover, the positive Bernoulli parameter 
indicates that ADAFs are prone to outflows although there have not been any 
self-consistent inflow/outflow solutions (Narayan Sz Yi 1995a). The vertically 
integrated equations and the self-similar solutions do not introduce any serious 
errors in flow dynamics since the height integration is a good approximation 
(Narayan & Yi 1995a). 
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2.3. Cooling and Heating Mechanisms 

In ADAFs, due to low radiative cooling efficiency of electrons and rather inef- 
fective ion-electron coupling, which is taken to be the Coulomb coupling, ions 
are nearly virialized (Narayan & Yi 1995b) 

Ti ~ 2 x \Q l2 f3r- l K. (11) 

Electrons' energy balance is maintained as (Mahadevan & Quataert 1997) 

ds de , m dn „ , . . 

^m = ' m m- Kr m = q +Sq " 9 (12> 

where kT e v(dn/dR) = q com P ress i s the compressive heating and (5(7 + is the direct 
viscous heating on electrons with 5 ~ m e /m p ~ 10 -3 (Nakamura et al. 1996). 
qie > qCompress £ Qr ^ > qj q 2 anc } p V (de/dR) m q %e — q~ describes the electron 
energy balance, cf e < q com P ress occurs when m ^ 10~ 4 a 2 and the energy balance 
for electrons becomes pv{de/dR) ~ q com P ress ; which is appropriate only for very 
low luminosity systems. < q com P ress [ s realized only for <5 ^ 1CT 2 which 

implies that the direct viscous heating is uninteresting in most practical cases 
(e.g. Mahadevan & Quataert 1997). 

For convenience, we introduce some physical scalings; mass m = M/Mq, 
radius r = R/R s (R s = 2GM/c 2 = 2.95 x 10 5 m cm), accretion rate m = 
M/M Edd (M Edd = L Edd /0.1c 2 = 1.39 x 10 18 m g/s where L Edd is the Eddington 
luminosity). The equipartition mag netic field B 2 /8ir = (1 - /3)pc 2 with (3 = 0.5 
and / = 1 where (3 is the ratio of magnetic to total pressure. The self-similar 
solution gives the physical quantities in terms of the dimensionless variables 
defined here. 

vm-lx W 10 ar~ 1/2 cm/s, (13) 
Q « 3 x 10 4 m-V~ 3/2 s _1 , (14) 
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1 x lO^r" 1 cm 2 /s 2 , (15) 

n e « 6 x 10 19 a -1 m -1 mr- 3/2 cm- 3 , (16) 

B^8x 10 8 a- 1 / 2 m- 1 / 2 r- 5 / 4 m 1 / 2 G, (17) 

r es « 24a" 1 r" 1/2 m, (18) 

9 + w 5 x 10 21 77?T 2 r~ 4 m erg/s/cm 3 , (19) 

where n e is the electron number density and r es is the electron scattering depth. 

Various (electron) cooling processes give rise to distinct spectral components 
(Narayan & Yi 1995b). The integrated cooling rate Q ~ J dVq where q is 
the cooling rate per unit volume and J dV denotes integration over the entire 
accretion flow. The total electron cooling rate is 

Q e — Qsync ~\~ Qcompt ~i~ Qbrem (^0) 

where Q syn c is the synchrotron cooling rate which gives rise to spectral emission 
components in radio, IR, or optical/UV depending on the mass m and accretion 
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rate rh. Qcompt 18 the Compton cooling rate which is mainly responsible for op- 
tical/UV/soft X-ray emission. Qbrem is the bremsstrahlung cooling contributing 
to X-ray and soft gamma-ray emission. If ADAFs in the inner regions around ac- 
creting compact objects are surrounded by the optically thick disks, optical/UV 
emission from cool disks is expected. Energetic protons in ADAFs may result in 
high energy gamma-rays (Mahadevan et al. 1997). Similar radiation processes 
in zero angular momentum spherical accretion have been extensively discussed 
(e.g. Melia 1992 and references). 

2.4. Critical Quantities 

ADAFs exist when accretion rates fall below a certain critical rate M cr u = 
rn cr itMEdd- Such a critical rate arises because there exists a maximum accretion 
rate above which heating could be balanced by radiative cooling without any 
necessity of advective cooling (Rees et al. 1982, Narayan & Yi 1995b, Abramow- 
icz et al. 1995, Narayan et al. 1998). In the case of the single temperature 
case, i.e. T e = T; L cx r -1 , assuming an optically thin, bremsstrahlung cooling 

(good for r > 10 3 ), we get q + cx rrT 2 r~ A rh cx rh and q~ = q^ rem oc p 2 T^ 2 ~ 
p 2 T cx a~ 2 m~ 2 r~ 7 / 2 rh 2 cx rh 2 or from q + ~ q~ we get the critical accretion rate 
rh cr it oc a 2 r~ 1 / 2 . In the case of the single temperature with T e = Tj / T V { r , 
assuming synchrotron and Compton cooling, the critical accretion rate becomes 
rhcrit ~ 10~ 4 a 2 . These well motivated critical rates are too low to be of practical 
interest (Esin et al. 1997). Esin et al. (1997) found that the critical rate de- 
duced from observed spectral transition in soft X-ray transients is much higher 
than the above rates. 

In the two-temperature ADAFs, the bottleneck in energy transfer from ions 
to electrons define another critical rate which is good for r ^ 10 3 (e.g. Narayan et 
al. 1998). That is, using q + cx rh and q~ cx q ie cx rh 2 and equating the two rates, 
q + = q %e , gives a critical accretion rate rh cr n ~ 1 x 10 3 (1 — /)e'a 2 /3~ 1 / 2 r~ 1 / 2 . 
Alternatively, tj e = t acc R/v gives rh cr it ~ 0.3a 2 . In sum, rh cr it ~ a 2 for 
r ^ 10 3 and rh cr u oc a 2 r~ l l 2 for r ^ 10 3 , which depicts the actual radial 
dependence of the critical accretion rate. 

It is interesting to point out that there exists a critical a (Chen et al. 1995). 
For a < a cr u ~ r, rh cr i t exists while for a > a cr u ~ r, rh cr i t doesn't exist. 

2.5. Some Recent Works on Heating Ions 

Bisnovatyi-Kogan & Lovelace (1997) recently claimed that large electric fields 
parallel to magnetic fields can accelerate electrons and hence bypassing the bot- 
tleneck in energy transfer from ions to electrons, which could rule out ADAFs as 
a possible accretion flow type. However, such a possibility is realized only when a 
substantial resistivity on microscopic scale exists. This requires a small magnetic 
Reynolds number. Their proposal to use the macroscopic turbulent resistivity 
is not applicable on microscopic scales as pointed out by Blackman (1998) and 
Quataert (1998). Blackman (1998) argued that the Fermi acceleration by large 
scale magnetic fluctuations associated with MHD turbulence may lead to pref- 
erential ion heating and hence two-temperature plasma. This heating is not 
applicable to non-compressive Alfvenic turbulence which is most likely to be 
more important than the compressive mode. For weak magnetic fields substan- 
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tially weaker than equipartition fields, Alfvenic component of MHD turbulence 
is dissipated on scales of proton Larmor radii (Gruzinov 1998, Quataert 1998). 
This mechanism favors ion heating and two-temperature plasma. For strong 
fields (i.e. near equipartition), the Alfvenic turbulence cascades to scales much 
smaller than proton Larmor radii and can directly heat electrons, which could 
cast doubt on ADAFs with equipartition strength magnetic fields. That is, there 
is a possibility that equipartition plasma doesn't allow two-temperature plasma. 
This issue may ultimately be resolved by observed spectra. 

2.6. ADAF Luminosity 

In ADAFs, the observed radiative luminosity Ladaf = / L{E)dE ~ J q ie dV or 
in terms of the ADAF radiative efficiency tjadaf, Ladaf = ijadafMc 2 where 
VADAF = f]eff x 0.2ma~ 2 ocmoc M/M (Narayan & Yi 1995b). In contrast, the 
thin disk luminosity L t h in disk ~ VeffMc 2 with r] e ff ~ 0.1 (e.g. Frank et al. 
1992). 

3. Some Physical Issues 
3.1. Global Solutions 

The self-similar solution considered so far applies to regions far from the inner 
and outer boundaries where physical scales of the system demand the dynamical 
equations to deviate from self-similarity. A comprehensive summary of the New- 
tonian case is found in Kato et al. (1998). The task of finding the global solution 
is to find an eigenvalue j (specific angular momentum accreted by the central ob- 
ject) with proper boundary conditions. The proper boundary conditions are (i) 
outer thin disk matching inner ADAF, (ii) sonic point for the transonic ADAF, 
and (iii) vanishing torque at the inner boundary. Solving the radial momentum 
equation, angular momentum equation, and energy equation, along with the 
implicit continuity equation for v,£l,c s ,p and eigenvalue j, global solutions are 
found. The self-similar solution is a good approximation for a wide range of 
radii between the inner and outer boundaries. The pseudo-Newtonian potential 
case has been solved by Matsumoto et al. (1985), Narayan et al. (1997), Chen 
et al. (1997). The major findings are as follows, (i) For a ^ 0.01, inefficient 
angular momentum transport results in slow radial accretion and a wide radial 
zone of super-Keplerian rotation. In the super-Keplerian rotation region, a thick 
torus-like structure with funnel around the rotation axis forms and the pressure 
profile shows a maximum, which is reminiscent of the ion torus model (e.g. Rees 
et al. 1982). (ii) For a ^ 0.01, efficient angular momentum transport and rapid 
accretion v oc avx occur. There exist no pressure maximum and the accretion 
flows are quasi-spherical. 

In the relativistic case with a spinning black hole (Abramowicz et al. 1996, 
Peitz & Appl 1997, Gammie & Popham 1998, Popham & Gammie 1998), the 
Newtonian and pseudo-Newtonian results are largely confirmed for r ^ 10. For 
r ^ 10, however, significant differences and spin effects are seen. Detailed calcu- 
lations of emission spectra have been carried out by Jaroszynski & Kurpiewski 
(1997). 
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It has been an issue whether shocks form in transonic accretion flows. In the 
steady calculations, shocks are not seen and in the time-dependent calculation of 
Igumenshchev et al. (1996), no shocks have been seen. Manmoto et al. (1996)'s 
time-dependent calculation shows shock-like steepening in density waves, which 
could be responsible for ADAF variabilities. 

3.2. Stability 

The geometrically thin disks could be thermally and viscously unstable under 
certain circumstances (e.g. Frank et al. 1992). ADAFs are primarily stable 
agains thermal and viscous perturbations: (i) ADAFs are stable against long 
wavelength perturbations (Narayan & Yi 1995b, Abramowicz et al. 1995). (ii) 
ADAFs are marginally stable against short wavelength perturbations in the 
single temperature case (Kato et al. 1996, 1997, Wu 1997). (iii) ADAFs are 
stable against short wavelength perturbations both thermally and viscously in 
the two temperature case (Wu & Li 1996). The stability of ADAFs have led to 
an argument that the thermally unstable Shapiro-Lightman-Eardley disk (Piran 
1978) is a spatially transitional accretion flow linking the thermally unstable 
outer thin disk to the stable inner ADAF. In the time-dependent calculations, 
small wavelength perturbations in the single temperature ADAFs have been 
observed to grow while moving inward. The growth rate is however not rapid 
enough to affect the steady global structure (Manmoto et al. 1996). 

4. X-ray Transients 

4.1. Black Hole Systems 

Detailed spectral fitting has been tried for black hole systems such as A0620-00, 
V404 Cyg (Narayan et al. 1996), Nova Mus 1991, Cyg X-l, GRO J0422+32, 
GRO J1719-24 (Esin et al. 1997), and 1E1740.7-2942 (Vilhu et al. 1997). The 
main physical parameters used in the spectral fitting are M, m, a, (3, and rt r , 
where the last one is the radius of accretion flow transition from outer thin 
disk to inner ADAF. The spectral fitting assumes that (i) the outer thin disk 
is joined to the inner ADAF and (ii) the outer thin disk is unstable against 
disk instability. The instability causes the heating/cooling waves to propagate 
inward, resulting in delays between different emission components (Lasota et al. 
1996, Hameury et al. 1997). The transition radius r% r is crucial in determining 
spectra (Lasota et al. 1996, Narayan et al. 1996,1998, Esin et al. 1998). 

4.2. Thin Disk - ADAF Transition 

ADAFs exist when q~ < q + which in principle determines rt r . In the single 
temperature, bremsstrahlung dominated case, 

q* is oc mr 2 r~ A m oc r -4 (21) 

Ibrem <* P^ 1 ' 2 K a~ 2 m~ 2 r' 1 1 2 m 2 oc r~' 7/2 (22) 
and q + = q~ gives (Honma 1996) 



r, r ~3x 10 2 (a 4 /m 2 ). 
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(23) 



In more realistic cases, details of transition are unknown. For instance, Honma 
(1996) considers the radial turbulent heat diffusion near the interface between 
the thin disk and the ADAF. Spectral fitting gives a different result on r tr (Esin 
ct al. 1997). The issue of the accretion flow transition still remains unresolved. 

4.3. Neutron Star Systems 

Neutron star transients have not been well studied using the ADAF models 
mainly due to the lack of sufficiently well-developed physical understanding of 
accretion flows near the neutron star surface. The radiation feedback from 
the soft radiation emitted by the stellar surface makes the spectral calculations 
considerably more complicated (Narayan & Yi 1995b, Yi et al. 1996). Moreover, 
it is unclear whether the spectral transition similar to that of black hole systems 
exists in neutron star systems. 

4.4. Accretion-Powered X-ray Pulsars 

Spectral transition in neutron star systems could be potentially very difficult to 
detect if thermalization near the neutron star surface occurs rapidly (e.g. Yi et 
al. 1996). Accretion-powered X-ray pulsars could however provide an observable 
signature of the temporal accretion flow transition. In accretion-powered pulsars 
such as 4U 1626-67, GX 1+4, and OAO 1657-415, abrupt torque reversals have 
been observed, which are extremely difficult to explain within the conventional 
models involving disk-magnetospheric interaction (Chakrabarty et al. 1997). 
The difficulties arise mainly due to (i) short reversal time scales, (ii) nearly 
identical spin-up and spin-down rates, (iii) small X-ray luminosity changes, (iv) 
significant spectral transition reported in 4U 1626-67 (Vaughan & Kitamoto 
1997), and (v) torque-flux correlation observed in GX 1+4 (Chakrabarty et al. 
1997). 

The observed phenomenon is well accounted for by the accretion disk tran- 
sition triggered by a gradual, small amplitude modulation of mass accretion 
rate. The sudden reversals occur at a rate ~ W 16 ~ 17 g/s when the accretion flow 
makes a transition from (to) a primarily Keplerian flow to (from) a substantially 
sub-Keplerian, radially advective flow (Yi et al. 1997, Yi & Wheeler 1998). The 
proposed transition model naturally shows that (i) the transition time scale is 
likely to be shorter than days, (ii) the required accretion rate change is at the 
level of a few xlO percent, and (iii) the abrupt reversal is a signature of a pulsar 
system near spin equilibrium with small mass accretion rate modulations near 
the critical accretion rate on the time scale of years. Other possible explanations 
for the spectral transition and the torque-flux correlation have been suggested 
(Nelson et al. 1997, van Kerkwijk et al. 1998) with varying difficulties in ex- 
plaining the above observational facts. The accretion flow transition is similar 
to those in black hole transients and cataclysmic variables, which strongly sug- 
gests a common physical origin (Yi et al. 1997). The transition time scale is 
tthermai ~ (a^x) -1 or t vis ~ Rjoic s ~ (R/ H)t t hermai , where the latter time 
scale becomes ~ 10 3 s for a ~ 0.3, R ~ 10 9 cm, M ~ 10 16 <7/s. This time scale is 
short enough to induce transition on a time scale of a day or less. 

If the accretion flow's thermal pressure is a substantial fraction £ of the 
dynamical pressure, i.e. c 2 s ~ ^R 2 VL 2 K , the thickness of the flow H ~ £ 1 / 2 i?, 
the radial accretion speed vr ~ —a^RQx, and the angular rotational frequency 
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fi ~ (1 - 5£/2 - a 2 ^ 2 /2) 1 / 2 ^^ = AQ K where £ -> is the Keplerian limit and 
A < 1 for ADAFs and ^4 = 1 for thin Keplerian disks. When A < 1 occurs, 
R' c = A 2 / 3 R C , where R c is the Keplerian corotation radius and R' c is the sub- 
Keplerian corotation radius. The accretion flow is truncated at a radius R' a and 
using TVq = AM (GM*R' ) 1 1 2 the torque on the star becomes 

= 7 1- (8/7) {R'jR' c fl 2 
% 6 l-{R'jR' c fl 2 1 ) 

which pushes the pulsar's spin toward an equilibrium spin P' eq = P eq j 'A > P eq 
where / denote quantities after transition from thin disk to ADAF. In this picture, 
torque reversal is expected if 

B.~5x 10 n Ll%ptf G (25) 

where L x ^q = L x /10 erg/s is the X-ray luminosity and -P*,io = P*/Ws is the 
pulsar spin period. Observed quasi-periodic oscillation (QPO) periods tightly 
constrain the proposed model. Yi & Grindlay (1998) discuss some possible 
implications on spin-up of LMXBs to MSPs when the accretion flows are ADAFs 
in these systems. 

4.5. Energetic Protons: Lithium Production 

Energetic ions present in ADAFs are capable of nuclear spallation. Lithium 
production in ADAFs in binary systems and ion tori in galactic nuclei has been 
discussed by Ramadurai & Rees (1985), Jin (1990), Yi k Narayan (1997). Using 
the self-similar solution, the relevant physical quantities are as follows: number 
density of protons nu ~ 6 x 10 20 m _1 mr _3 / 2 cm -3 , number density of a particles 
n a ~ 5 x 10 19 m~ 1 mr -3 / 2 cm -3 , energy per nucleon E ~ 300r~ 1 MeV, and the 
radial accretion speed vr ~ 2 x 10 9 r -1 / 2 cm/s. The production of 7 Li domi- 
nates and the production cross-section a+(E) ~ 100(E / 100 MeV)~ 2 mbarn for 
E > 8.5MeV. Continuous production of 7 Li within the accretion flow leads to 
increase in Lithium abundance according to 

An Li 1 n 2 

= oC+l-^K — At f i ow (26) 

where Atfi ow = AR/vr. The enrichment results in the terminal abundance 

nu/nu = J d{nLi/nu) ~ 0.1m or the total Lithium mass production rate Mli = 
2 x lO" 8 mm 2 M /V. 

It is expected that 7 Li enrichment occurs around NSs and BHs containing 
ADAFs, which has been seen in recent precision spectroscopic measurements of 
7 Li in V404 Cyg, A0620-00, GS 2000+25, Nova Mus 1991, Cen X-4 (Martin et 
al. 1994 and references therein). In contrast, WD systems show no such effect, 
which indicates that only at r ~ 1 relativistic energies are reached while r ~ 10 3 
is the inner most radius in the WD systems, which gives E <C lMeV. 

4.6. Thermalization of Particles 

Although we have adopted the thermal temperatures for ions and electrons, it 
is not proven that particle energy distributions are adequately approximated by 
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thermal distributions. Recent investigations (e.g. Mahadevan &; Quataert 1997, 
Quataert 1998, Blackman 1998) have suggested some limited information on this 
unresolved issue. 

In most of the ADAF models, protons and ions are energized by viscous 
heating which is mostly unspecified. Alfvenic turbulence (which does not re- 
sult in strongly non-thermal distributions) and Fermi acceleration (leading to 
power-law tails) have been considered. Coulomb collisions and synchrotron ab- 
sorption do not lead to rapid thermalization of protons. So far it appears that the 
acceleration mechanism itself determines the proton energy distributions. Ther- 
malization of electrons could occur more easily (Ghisellini & Svensson 1991). 
Coulomb collisions can thermalize electrons for m J> 10 _2 a 2 . Synchrotron self- 
absorption leads to thermalization when rh 10 _5 a 2 r. The electron energy 
distributions directly affect radio emission spectra. 

Recently Mahadevan (1998) pointed out that low frequency v ^ 10 9 Hz ra- 
dio spectrum of Sgr A* could be contributed by electrons and positrons produced 
by charged pions from proton-proton collisions and that neutral pion production 
could account for tentative detection of gamma-rays in the direction of Sgr A*. 
In both cases, particle distributions need to be strongly nonthermal and the 
results depend very sensitively on high energy tails. 

5. Galactic Nuclei 

ADAFs may exist in galactic nuclei including the Galactic center source Sgr 
A*. For our discussions, we define mj = m/10 7 , m_3 = m/10~ 3 , and R s = 
2GM/c 2 = 3 x 10 12 m 7 cm. 

When the accretion flow is a thin disk, the luminosity L = rjMc 2 with 
the efficiency r\ = r] e ff ~ 0.1 (e.g. Frank et al. 1992). Although the total 
luminosity is high, the emission temperature or the disk temperature ~ 
6 x 10 6 77i 7 1 ^m'J_^r~^^K is too low to account for X-ray emission. The disk 
luminosity Ldisk ~ 1 x lQ 42 rnjrn^erg/ s is expected to occur in optical/UV/soft 
X-ray. X-ray emission and radio emission are usually accounted for by optically 
thin corona and radio jets. For the latter, the estimated radio power is Lj e i ^ 
1 x 10 42 a 2 r}j et rn7m-3erg/s where a is the black hole spin parameter and f]j et is 
the jet radiative efficiency. 

In ADAFs, relevant physical quantities scaled for galactic nuclei are equipar- 
tition magnetic field B ~ 1 x 10 4 m 7 1 rh_ 3 2 r~ 5 / 4 G, electron scattering depth 
T es ~ 5 x 10 -2 m_3, the ion temperature Tj ~ 2 x 10 12 r~ 1 K, and the electron 
temperature T e ~ 5 x 10 9 K. ADAFs are expected to exist when the mass 
accretion rate falls below fn cr it = M cr i t / MEdd ~ 0.3a 2 ~ 10~ 3 — 10~ 2 . Radio 
emission is easily explained by the synchrotron emission with the characteristic 
synchrotron emission frequency (Yi & Boughn 1998ab) 

v sync ~ 1 x 10 12 mf 1/2 my 3 V 5 / 4 T 2 9 Hz (27) 

where T e g = T e /W 9 K ~ 5. The highest synchrotron radio emission frequency is 
Vmax = v sync (r ~ 1) ~ 3 x 10 13 m 7 1 ^ 2 m 1 /^ Hz which comes from the inner most 
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region of ADAF near the black hole horizon. The radio luminosity 

L R ~ uLfT ~ 2 x l^x%\TT m ^m^ erg/s (28) 

where xm3 = ^m/10 3 is a dimensionless synchrotron self-absorption parameter 
and vio = v/W 10 Hz. ADAF radio emission could be tested by the radio source 
size - frequency relation. ADAFs predict such a relation 

0{y) ~ 2m 3 7 /5 m 2 ^u w 4/5 (D/WMpcy 1 fias (29) 

where the angular size 9{v) becomes ~ mas for distance scales D ~ Wkpc. 

Optical/UV/X-ray in ADAFs arise from inverse Compton scattering of ra- 
dio synchrotron photons and hard X-rays are from bremsstrahlung and mul- 
tiple Compton scattering. For m ^ 10~ 3 , X-ray emission is dominated by 
the bremsstrahlung emission, L x ~ L x rem oc mm 2 and the radio luminosity 

Lr oc m 8 / 5 m 6 / 5 , which suggests Lr oc mL^/ 5 . For m ^ 1CT 3 , X-ray emission 
gives the luminosity L x ompt oc fn 7 ^ +N with N > 2. In this case, we expect 

Lr oc mLi^ +1 '. Yi & Boughn (1998ab) have derived and tested the radio/X- 
ray luminosity relation for ADAFs using L x = L x (2 — WkeV); 

L R ~ 10 36 m 7 (^/15GFz) 7 / 5 (L :c /10 40 er 5 /s) 2: er 5 /s (30) 

where x ~ 1/5 for m ^ 10~ 3 and x ~ 1/10 for m > 10~ 3 or L R ^ adv / 'L x ^ v oc 
mL xMv 

ADAFs are likely to drive jets/outflows (Narayan & Yi 1995a). If jets are 
powered by black hole's spin energy (e.g. Frank et al. 1992), 

LR,jet/L Rja dv ~ 4 x 10 5 a 2 7] jet my 1/5 m 1/ 3 5 (31) 

and L Rtjet /L Xjadv oc a 2 mL^ adv are expected. That is, L Rtjet > L R , adv for a > 
2 x 10 -3 77~y 2 . It is often argued that radio-loud nuclei have a ^ 1. If a galactic 

nucleus contains a thin disk and a jet L R j et / 'L x ^isk ~ a?£jet/Veff ~ 0(1) is 
expected. 

Characteristic ADAF emission spectra are determined primarily by m and 
weakly affected by the black hole mass M. Any combinations among L x , Lr 
and M give useful information on the nature of emission from galactic nuclei. 

5.1. Galactic Center Source Sgr A* 

Galactic center radio source Sgr A* is a prime candidate for an ADAF (Rees 
1982, Narayan et al. 1995, 1998, Manmoto et al. 1997). Sgr A*, which is at 
the dynamical center of the Galaxy, appears to contain a massive black hole 
with mass ~ (2.5 ± 0.4) x 10 6 M Q . Wind accretion from nearby IRS 16 wind 
is expected to provide a mass accretion rate M }t a few xW~ 6 MQ/yr. With 
the conventional ~ 10% efficiency, such a high accretion rate would correspond 
to a luminosity of O.lMc 2 ^ 10 40 erg/s which is some 3 to 4 orders of magni- 
tude larger than the observed radio to gamma-ray luminosity of ^ 10 37 erg/s. 
Moreover, a standard thin disk would give peak emission in near infrared but 
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Menten et al. (1997)'s 2.2 micron upper limit rules out this possibility. These 
facts strongly suggest that an ADAF is present in Sgr A*. 

Spectral fitting based on the observed M and the estimated m adequately 
accounts for the observed emission seen from radio to hard X-ray. Bower & 
Backer (1998) measured the intrinsic source size ^ 0.48 mas at 7mm, which 
at 8.5kpc gives the linear size of 4.1AU and the lower limit on the brightness 
temperature 4.9 x 10 9 K. Both the size and temperature measurements are con- 
sistent with the ADAF models. Inverted radio spectrum with sharp cutoff along 
with no jet like elongations near Sgr A* are also consistent with the ADAF pre- 
dictions. X-ray constraints are rather uncertain. ROSAT 0.8-2.5 keV luminosity 
is ~ 1.6 x 10 3A erg/s (Predehl & Trumper 1994) and ASCA 2-10 keV luminosity 
is < 4.8 x 10 35 erg/s (Koyama et al. 1996). Both constraints are easily satisfied 
by an ADAF. 

5.2. NGC 4258 

NGC 4258 almost certainly has a central black hole with mass M = 3.5 ± 
0.1 x 10 7 M Q , which is concentrated within 0.13pc (for a distance of 6.4Mpc) 
of the dynamical center. This source has an observed 2-10 keV X-ray lumi- 
nosity of L x = 4 x 10 40 erg/s (Makishima et al. 1994). Optical luminosity 
L op t ~ \Q A2 erg/s (Wilkes et al. 1994) provides an additional constraint. 22 
GHz continuum emission (after subtraction of jet component) has not been de- 
tected with a 3cr upper limit of 220/xJy (Herrnstein et al. 1998) or a luminosity 
upper limit Lr{22GHz) < 2.4 x 10 35 erg /s. The non-detection of the core at 
22 GHz could imply that m ~ 10~ 2 and r tr ~ 30 if an ADAF exists in NGC 
4258 (Gammie et al. 1998, cf. Lasota et al. 1996). Such a constraint is highly 
suspect due to a possibility of strong variabilities in radio emission from ADAFs 
(Blackman 1998, Ptak et al. 1998, Herrnstein et al. 1998). 

5.3. M60, M87, NGC1068, and M31 

Di Matteo & Fabian (1997b) has attempted to fit M60 emission spectra with 
an ADAF under the assumption that the accretion rate is close to the Bondi 
accretion rate with M ~ 10 9 Mq. Due to uncertainties and lack of flux measure- 
ments, a definitive conclusion as to whether an ADAF exists needs more data 
in X-ray or other wavebands. Reynolds et al. (1996) claimed that an ADAF 
model spectrum for a blackhole mass M = 3 x 10 9 M Q and the accretion rate 
m ~ 10~ 3 accounts for the observed fluxes of M87. However, such a conclusion 
is highly suspect because the extended radio emission component has not been 
properly removed. ADAFs themselves do not produce extended radio emission. 
NGC1068 is an obscured Seyfert with a very high obscuration-corrected X-ray 
luminosity X-rays are most likely to be dominated by scattering and the X-ray 
luminosity may be too bright for ADAFs based on the estimated black hole mass 
~ a fewxlO 7 M (Yi & Boughn 1998ab). Yi & Boughn (1998b) also considered 
M31 which has a central black hole of mass M = 3 x 10 7 M Q . The observed radio 
luminosity is too low for the observed X-ray luminosity if an ADAF is assumed 
around the black hole. It is highly likely that the X-ray luminosity is dominated 
by binary sources in the nucleus while the radio emission may be due to a very 
weak ADAF. 
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5.4. X-ray Bright Galactic Nuclei 

Yi Sz Boughn (1998ab) and Franceschini et al. (1998) applied the ADAF model 
to a small sample of X-ray bright galactic nuclei which have black hole mass es- 
timates. Since ADAFs are most relevant for low luminosity, hard X-ray sources, 
faint, hard X-ray galactic nuclei are likely hosts of ADAFs (Fabian & Rees 1995, 
Di Matteo & Fabian 1997a, Yi & Boughn 1998ab). Hard spectrum, faint X-ray 
sources could contribute significantly to the diffuse X-ray background. 50% of 
2-10 keV XRB could be accounted for by ADAF sources with the comoving 
density of 3 x 10~ 3 Mpc~ 3 for L x ~ W 41 ergjs, which is comparable to the local 
density of the L* galaxy (Di Matteo & Fabian 1997a, Yi & Boughn 1998a). 
However, unless M ^ 1O 9 M , L x ^ 10 40 erg/s would be already too bright 
for ADAFs to account for the observed X-ray background. This is because at 
high luminosities X-ray emission is dominated by the Compton scattering which 
result in X-ray spectra much different from the background spectrum similar 
to the bremsstrahlung-dominated X-ray spectra. That is, in relatively bright 
ADAFs, high m's correspond to the Compton-dominated cooling regime. Al- 
though a significant clumping of ADAF gas would enhance bremsstrahlung over 
Compton (Di Matteo &: Fabian 1997), such a possibility is difficult to realize. 

So far, we have argued that ADAF models in galactic nuclei are testable 
due to distinguishing characteristics of ADAFs. ADAFs with low radiative ef- 
ficiency and high temperature are likely for massive black holes accreting at 
accretion rates ^ l§~ 2 Mem- Hard X-ray emission and inverted spectrum radio 
emission from compact core are expected. There exists a characteristic radio/X- 
ray luminosity relation as shown by Yi & Boughn (1998ab). For known black 
hole masses, existence of hot ADAFs can be tested by radio/X-ray observations. 
Black hole masses could be estimated based on radio/X-ray luminosities. ADAF 
sources could however contain jets/outflows which can contribute to radio emis- 
sion. Depending on the level of radio activity and existence of extended radio 
emission features, galactic sources could be classified (Yi k, Bough 1998ab). Such 
a classification can be quantified in a manner similar to that adopted for Galac- 
tic X-ray sources. In fact, there exist interesting spectral similarities between 
Galactic binary X-ray sources and galactic nuclei. 

We define X-ray bright galactic nuclei (XBGN) as galactic nuclei with X- 
ray luminosities in the range 10 40 5; L x ^ 10 42 erg/s which is sub-luminous 
compared with the more powerful active galactic nuclei (AGN) which generally 
have L x ^ 10 43 erg/ s. Most of XBGN are expected to overlap with emission 
line galaxies with L x ~ 10 39 — 10 42 erg/s. However, some of the low luminosity 
Seyferts with L x ^ 10 42 ergjs cannot be ruled out. Based on our discussions of 
ADAFs, for L x ~ 10 41 ergjs, 

L fi ~ 4 x 10 36 (M SH /3 x 1O 7 M ) ergjs (32) 

at 20 GHz with the characteristic inverted radio spectrum I v oc u 2 ^ . These 
sources at distances ~ 10(M/3 x 10 7 M & ) l / 2 Mpc should be detected as ~ lmJy 
point-like radio sources (Yi & Boughn 1998a). If X-ray and radio are indeed 
from ADAFs, the black hole masses can be estimated (Yi & Boughn 1998b). 

Yi and Boughn (1998ab) proposed the source classification based on the 
known black hole masses and the ADAF radio/X-ray luminosity relation. Adopt- 
ing Sgr A*, NGC 4258, NGC 1068, NGC 1316, NGC 4261, and NGC 4594 as 
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fiducial sources, a statistically incomplete sample of XBGN are classified into 
radio-loud XBGN and radio-quiet XBGN. The former show that the observed 
radio luminosity L Rfibs ~ L RJet > L Radv where L RJet and L R adv are the ex- 
pected radio jet luminosity and ADAF radio luminosity, respectively. These 
sources are expected to show extended radio emission, unlikely to have strongly 
inverted radio spectra, and may have compact ADAF radio emission from com- 
pact cores separate from the extended emission components. The latter show 
that L R o t, s ~ L R ^ adv and that the dominating emission components are compact 
cores with inverted spectra. 

Kellermann et al. (1998) and Falcke (1998) show that jet-like radio emission 
features are common among AGN and emission line galaxies. Surprisingly, even 
in radio-quiet sources, elongated radio emission features are sometimes seen on 
small scales, which could imply that some type of jet/outflow activities are very 
common in galactic nuclei regardless of their large scale radio activities. In order 
to resolve this issue, high resolution radio measurements for nearby (< lOMpc) 
sources are crucial. Hard X-ray emission and inverted spectrum, compact radio 
emission are very likely to be found closely correlated. 

5.5. QSO Evolution 

Yi (1996) has suggested that transition of accretion flows from thin disks to 
ADAFs at rh cr it ~ 0.3a 2 could account for the observed sudden decline in the 
number of bright QSOs at redshift z ~ 2 and downward (see also Fabian & Rees 
1995). Once ADAFs set in, the luminosity evolves according to 

L = Ladaf ~ 30m x L E dd (33) 

where x ~ 2 and m = M /Mem oc M/M. the last expression implies that even 
when the mass accretion rate M is kept constant, m decreases merely due to 
the growth of black hole mass. 

For instance, in a flat universe with no cosmological constant, for the initial 
black hole mass M = Mi at z = z% and M = constant, when m < m cr i t and 
M/H > Mi 

L(z) oc (1 + z) 3 ^- 1 )/ 2 [(1 + z)- 3 / 2 - (1 + Zi )- 3 / 2 ] (34) 

L(z) oc {l + z) K ^ (35) 
3(x - 1) 



or 



with 



K(z) 



1 



1 



(36) 



(l+Z i )/(l+^)]3/2_l_ 

which shows that the luminosity declines with a power-law similar to that seen 
in observations (Yi 1996 and references therein). The epoch at which ADAFs 
set in for m = rhi ~ 1 is (i.e. m = m cr i t first occurs) 



l + z c 



(¥) ( J --^)+( i+ ^ 3/2 

V t J \m cr it oj 



-2/3 



(37) 



where 5 = Mt Edd /Mi and t Edd = (Me^/M)" 1 = 4.5 x 10 7 (r/ e// /0.1)yr. The 
observed sudden decline of QSOs at z ~ 2 is naturally accounted for. 
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5.6. Outflows 



ADAFs are prone to outflows or jets (Narayan & Yi 1994,1995a) although a self- 
consistent inflow/outflow solution has not been found yet (cf. Xu & Chen 1997). 
It remains to be seen if a self-consistent inflow/outflow solution can account for 
compact and extended jet-like emission components in XBGN. 



6. Some Unresolved Issues 



ADAFs may exist in sources spanning many decades of masses of compact ac- 
creting sources. Some of the old outstanding issues, which are mostly concerning 
with low-luminosities and hard X-ray emission, are plausibly resolved by various 
versions of ADAF models. There exist however a number of unsolved problems 
in the ADAF framework. 

(i) The physics of accretion flow transition, temporal and spatial, remains 
unclear. The spatial transition (i.e. from outer thin disk to inner ADAF) is 
better understood than the temporal transition to a certain extent. However, 
there has not been an adequate explanation for r% r . It remains unsolved why the 
disk flow makes a transition to ADAF with little, if any, change of accretion rate, 
(ii) Even luminous systems show energetic X-ray emission which is absent in the 
thin disk models. It is often assumed that these sources have X-ray emitting 
coronae for which little physics is known. It is crucial to link the ADAF models 
to corona models with proper physical understanding, (hi) ADAFs emission 
could be highly variable. The issue of steady vs. non-steady ADAFs is directly 
related to the observed variabilities in ADAF candidate sources which show 
occasional non-detections, (iv) A number of plasma astrophysical issues remain 
to be solved. Is the two-temperature flow physically allowed? Are particles 
rapidly thermalized? What is the correct strength of magnetic helds responsible 
for synchrotron emission? (v) Observationally, faint X-ray sources, which could 
be seen by high resolution, high sensitivity experiments such as AXAF, should be 
studied in great details. Ultimately, ADAF-related issues and the very question 
on the future of ADAFs' are likely to be answered by observations. 
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